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Determining CP violation angle 7 with B decays into a scalar/tensor meson 
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We propose a new way for determining the CP violation angle 7 without any hadronic un- 
certainty. The suggested method is to use the two triangles formed by the decay amplitudes of 
B ± ->■ (D°,D ,D a C p)K*f 2) {l<i'S0). The advantages are that large CP asymmetries are expected in 
these processes and only singly Cabibbo-suppressed D decay modes are involved. Measurements 
of the branching fractions of the neutral Bd decays into DKq( 2 ) (1430) and the time-dependent CP 
asymmetries in B s ->• [D°,D°)M (M = / (980), / (1370), /a(1525), /i(1285), /i(1420), /ii(1180)) 
provide an alternative way to extract the angle 7, which will increase the statistical significance. 
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CP violation in the standard model (SM) originates 
from a single, irreducible phase in the 3x3 quark mixing 
matrix called the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. Precision test of its unitarity allows us to ex- 
plore the SM description of the CP violation and reveal 
any physics beyond the SM. One of the foremost tasks 
during the past decades has been to study the so-called 
(bd) unitarity triangle, the graphical representation of 
the condition stemming from the unitarity of the CKM 
matrix: V ud V* b + V cd V* b + V td V* b = 0. Its sides can 
be measured by leptonic and semileptonic meson decays, 
while the determinations of the angles (a, /3, 7), satisfying 
a + j3 + 7 = 180°, rely mostly on nonleptonic B decays. 

Our knowledge of the angle /3 to a large extent benefits 
from the gold-plated channel B — > J/ipKs and the cur- 
rent results already have a precision better than 1° [H,[2j]. 
The accuracy on the angle a is around 4°, thanks to 
the measurements of charmless tree dominated processes 
B -» tt + tt~,B — > p ± n T ,B -» p + p~ and B — > a^ir T . In 
contrast, results for the angle 7 are less accurate, with a 
precision of roughly 10°, which is one of the main sources 
of the current uncertainties on the apex of the unitary 
triangle. 

Since the angle 7 = arg(-V ud V* b / (V cd V* b )) is the rel- 
ative weak phase involving the decays b — > cus and 
b — » ucs, several methods on the basis of the decays 
B ± -> DK ± , with D being any admixture of D° and 
D°, have been proposed (for a review, see Ref. [H). 
The most productive ones are the Gronau-London-Wyler 
(GLW) method [HI], with D decaying into the CP eigen- 
states including tv°Ks, ir + ir~Ks, K + K~it + it ~, the 
Atwood-Dunietz-Soni (ADS) method 0, Q, using the 
Cabibbo-favored and doubly Cabibbo-suppressed D de- 
cay modes, and the Giri-Grossman-Soffer-Zupan (GGSZ) 
method [§] , which makes use of a Dalitz-plot distribution 
of the products of the multi-body D decays. All three 
methods are theoretically clean and do not require any 
time-dependent measurement. 

In the GLW method, the sensitivity of the CP asymme- 



tries to 7 is proportional to the ratio of the two interfering 
amplitudes, which is of the order 10%. The ADS method 
demands a detailed knowledge of the doubly Cabibbo- 
suppressed D decays, while the GGSZ method requires a 
Dalitz-plot analysis of multibody D decays. In this work, 
we propose a new method which is based on B — > DM 
decays with M being a light scalar/tensor meson. The 
proposed method has both advantages, namely on the 
one hand the interference and the CP violation in the cho- 
sen decay modes are sizable and on the other hand nei- 
ther doubly Cabibbo-suppressed D decays nor the Dalitz 
plot are needed. Among the various B decays into a 
p-wave scalar/tensor meson to be discussed, of partic- 
ular interest are the B ± -> (D°, D°, D% p )K*^ (1430) 
modes, where K*^(1A30) is a scalar (tensor) meson 
with J p = + (2 + ). The small (zero) decay constant of 
K$ (1430) (if| (1430)) compensates the large Wilson coef- 
ficient in the color-allowed amplitude, resulting in similar 
sizes for the decay amplitudes of B ± -> D°K*± 2) (1430) 

and B^ 1 —} D°K*^ (1430). As a consequence, there are 
large CP asymmetries. Measurements of branching ra- 
tios (BRs) of the neutral B d decays into D^T*^ (1430) 
and time-dependent CP asymmetries in B s —> DAI (M = 
/ (980), / (1370), ^(1525), /i(1285), /i(1420), %(1180)) 
provide an alternative way to extract the angle 7. For 
the sake of brevity, hereafter we use Kq 2 and /o, to ab- 
breviate i<:£ )2 (1430) and /o(980),/£(1525), respectively. 

All three methods [Hi] to extract 7 based on B — > 
(D° , D° , D(~, p )K ± use the information that the six de- 
cay amplitudes form two triangles in the complex plane, 
graphically representing the following identities 
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where the convention CP\D ) = \D°) has been adopted 
and D , (D°_ ) denotes the CP even (odd) eigenstate. The 
corresponding Feynman diagrams for these processes are 
given in Fig. [TJ Measurements of the decay rates of the 
six processes completely determine the sides and apexes 
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pothesis 
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A(B~ 



D°K, 
D°K, 
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-v ub v*c, 

-V cb V: s (C-T), (2) 



0) 



where C = G F f D a 2 (m 2 B - m^O^f K {m 2 D )/y/2, T = 



FIG. 1: Feynman diagrams for the color-suppressed contri- G F f K . ai (m 2 - m 2 )F^ D (m 2 ,)/V2, and G 



Z?X ( - 2) (1430) (a), B 



butions in the process B 

D°Ko(2) (1430) (b), and the color-allowed contributions in the 
B~ — > D° -^0(2) (1430) (c). In the diagrams (a, b), the specta- 
tor quark can also be a d or s quark, in which the light hadron 
consists of #5(1430), / (980) and /a(1525). 



of the two triangles, in particular the relative phase be- 
tween A(B~ -> D°K~) and A(B+ -> D°K+) is 2 7 . 

The shape of the two triangles is governed by two quan- 
tities 



rp is 



the 



'■' = \A(B- -> D°KJ)/A(B~ 
g [e l ~<A(B- -> D°KJ)/A{B- 



D°KJ)\, 



with Kj = K,K* 2 . The B~ -> L>°A- is both Cabibbo- 
suppressed and color suppressed. Thus the ratio ~ 
|KbK*s/(^ r cbKTs) a 2/ai| ~ 0.1 is small and in fact the 
world averages for the parameters Q 



= 0.107 ±0.010, 5f = (112 



5+12 
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indicate that the two triangles are squashed. Physical 
observables to be experimentally measured, defined as 
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have a mild sensitivity to the angle 7, and their values 
are expected to be Rqp± ~ 1 anci ^cp± ~ 0. 

Since the -?QJ(2) have the same flavor structure as the 
K meson, the relations given in Eq. (JlJ also apply to 
B ± -> (L> ,L» ,L>^p)K* ( ^ ) (1430). We wish to point out 
that, because of the suppression of the color-allowed de- 
cay amplitudes, the low sensitivity problem is highly im- 
proved and in particular large CP asymmetries are ex- 
pected. Although the A^^-emission diagram, as de- 
picted in Fig.QJc), has a large Wilson coefficient a\ ~ 1, 
the emitted meson is generated by a local vector or axial- 
vector current (at the lowest order in a s ), whose matrix 
element between the QCD vacuum and the Aq(A|) state 
is small (identically zero). 

A crude and model-dependent estimate of the ampli- 
tudes can be made with the help of the factorization hy- 



Fermi constant. The decay constant, defined via 
(K*-(U30)\s^u\0) = f KS p^., 

vanishes in the SU(3) symmetry limit and may get a 
nonzero but small value due to the symmetry break- 
ing effects. The current experimental data on t — > 
A'q~(1430)Sv places an upper bound [1(3] 

\f K5 \ < 107MeV, 

which is not very st ring ent. Adopting an estimate based 
on QCD sum rules [11[ 



-24MeV, or f K * = 36MeV, 



which contains a sign ambiguity, we find the relation 
2ai|/if*| ~ 0.2/0, with the D meson decay constant ex- 
tracted" from D~ -> fxD^. f D = (221 ± 18)MeV 
Using one set of results for the B — > AjJ form factors 
calculated in the perturbative QCD approach [l^] (cor- 
responding to fx* — 36 MeV), the B — » D form fac- 
tors from Ref. [l3| and 0,2 = 0.2, a\ = 1 we estimate 
C/T~ 1.2 and 



rf = \cv ub v;j[v cb v; s (c - t)]| . 

The corresponding BRs are roughly 

B(B- -> D°K*-) - 4 x 10~ 6 . 



(3) 



(4) 



Since the strong phase can not be computed at present, 
we take several ^benchmark values to illustrate the de- 
pendence of Rcp-l an d A c ° p+ in Fig. [2] In panels (a,b), 



k* 



2 is employed, and in panels (c,d) = 1. In 
the last two panels (e,f), we consider the case in which 

the ratio is not enhanced too much — 0.3. The 
solid (green), dashed (black), dotted (blue) and dot- 
dashed (orange) lines in diagrams (a,c,e) are obtained 

with 6g° = (30,60,120,150)° respectively, while the 
corresponding lines in diagrams (b,d,f) correspond to 

30, -60)°. The shadowed (light-green) 
region denotes the current bounds on 7 = 

(68112)° from 

a combined analysis of 5+ — > DK ± [2], in which the 
vertical (red) line corresponds to the central value. The 
CP odd quantities can be obtained similarly, for instance 

R-CP- — (-R-CP+) x K o 



S B ° = (30,60 



<S„° -^180° -S 
± 



Turning to the B ± — ¥ DK^ mode in which the ma- 
trix element of the vector and the axial- vector current be- 
tween the QCD vacuum and the A 2 state is zero, we find 
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FIG. 2: The dependence of Rq P+ and A^p + on 7. In 
panels (a,b), rg° = 2 is employed, in panels (c,d) rg° = 1 

and in panels (e,f) r^ — 0.3. The solid (green), dashed 
(black), dotted (blue) and dot-dashed (orange) lines in dia- 
grams (a,c,e) correspond to 5 B ° = (30,60,120,150)° respec- 
tively, while the corresponding lines in diagrams (b,d,f) cor- 
respond to 5g° = (30,60,-30,-60)°. The shadowed (light- 
green) region denotes the current bounds on 7 = (68j^5)° 
from a combined analysis of B — ¥ DK ± 3, in which the 
vertical (red) line corresponds to the central value. 



a vanishing color-allowed amplitude T. Accordingly, the 
ratio r B 2 is from the product of CKM matrix elements 
which is roughly 0.5. An estimate of the branching ratios 
can be made by using the data on the B — > J/^K 2 



B{B- -S- D°K* 2 
B(B 



J/^Kf) 



h 



V cb V* f,, N , 



0.8%,(5) 



with xk* being the ratio of the form factor products 
which is evaluated from a recent calculation of B — » K 2 
form factors (l4j ]: xk* — 0.5. The branching ratio 



B(B -> J^Kf) = 
from the data on B 

B(B~ -> 



(4.0 ± 2.4) x 10~ 4 




d°k;-) 



3 x 10" 



151 ] extracted 
gives 

(6) 



The method to use the two triangles formed by the 
six decay amplitudes for determining 7 is also valid in 
the neutral Bd decays into DKq° 2 , in which the tree am- 
plitude T is identically zero. The Kq 2 is self-tagging, 
thus no time-dependent measurement is required. Since 
the amplitudes involving D° and D arise from the same 
type of diagram, one expects that 5 B ° ~ 0. If true, 
the CP asymmetries A^p^ would be still close to but 
R C p± can largely deviate from 1. 



The long-distance contributions in the form of final 
state interactions (FSI) might change the factorization 
analysis in at least two aspects. First, FSI can give non- 
trivial strong phases to C and T which are zero in the 
factorization approach. Second, FSI might also modify 

the size of the amplitudes and the r B 0,2 . Despite these 
changes, no hadronic uncertainties will be introduced as 
the CKM matrix elements in the final state interactions 
are the same as the ones in Eq. @. To account for 
such effects, we also show in Fig. [2] the dependence of 
on 7 with different ratios of amplitudes: 



1 and 



K' 



0.3. The latter corresponds to the 



sign of Wilson coefficient a<i reversed namely ai 

K* 



-0.2. 



In this case, despite a small ratio r B ° = 0.3 the branching 



D Q K, 



q 2 ) can reach 10 



fractions B(B 

The Dq P meson in the final state can be recon- 
structed in the CP eigenstates, including the modes 
ir°Ks, 7r + 7T~ Ks, K + K~ , 7r + 7T~. These modes have quite 
large BRs, for instance, B(D° ->■ tt+tt-^s) - 3% [nj. 
The Kq 2 have significant decay rates into K tt, with 
B(K* 4 Kit) = (93 ± 10)% and B{K* Kir) = 
(49.9 ±1.2)%, and the final mesons are also easy to detect 
in experiments at hadron colliders. Moreover, since the 
CKM matrix elements for the Kq and K 2 are the same, 
no knowledge of the resonance structure in this method is 
required and therefore the angle 7 can be extracted with- 
out any hadronic uncertainty. Compared with the BR of 
B~ — > D Q K~, of order 10~ 6 , which is an unavoidable 
entry in the currently-adopted methods to determine 7, 
the summed BRs for the channels involving Kq and K 2 , 
of order 10 -5 , are comparable or even larger, and hence 
their measurements will not be statistically limited. The 
large amount of data accumulated by LHCb recently and 
in future will lead to a promising prospect of the proposed 
method. 

In the above discussion, we have neglected effects 
caused by the CP violation in D decays which is an- 
ticipated to be small in the standard model. Based 
on the 0.62 fb _1 of data collected in 2011, the LHCb 
collaboration [17| has measured the difference between 
CP asymmetries in singly Cabibbo-suppressed decays 
D° — > K + K~ and D° tt+tt-, AA cp = A CP (D° 
K+K-) - A C p(D° tt+tt-), given by 

AA CP = (-0.82 ± 0.21(stat.) ± 0.11(sys.))% , (7) 

where the first uncertainty is statistical and the second 
is systematic. Together with the CDF results [TH and 
previous world average from Heavy Flavor Averaging 
Group 19], the new world average for AAqp is found 
to be [2(| 



AA CP = -(0.645 ± 0.180)% 



(8) 



Although the new world-averaged AAqp is about 3.6ct 
away from zero, its magnitude is smaller than 1 percent. 
As a consequence, the CP violation effects in charm de- 
cays are less important in our method to determine 7, es- 
pecially when compared to large uncertainties in current 
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knowledge of 7. Moreover, since the direct CP violation 
in D° — > K + K~~ and D° — > 7r + 7r~ modes is expected 
to have opposite signs, part of the CP violation effects 
will cancel when both decay modes are used in the re- 
construction of D meson. 

Now we turn to the B s — >■ DM decays, whose Feynman 
diagrams are depicted in Fig. [T] with q = s. It is proposed 
in Ref. p| [2l[ that the time-dependent CP asymmetries 
in B s — > Dtp can be used to extract 7 and this method 
is ap plie d to a pure annihilation mode B s — > D^n^ in 
Ref. (Hi and modes like B s Dr)(r)') in Ref. [H|. In 
the example of B s — > Dfo, there are four decay modes 
having the amplitudes 

A(B S -> D°f ) = V ub V c * s A u A{B s -> D°/o) - V^V^, 
A(B S D°f ) = V cb V: s A 2 ,A(B s ->■ D°f ) = V c * b V us A 2 .(9) 

For each amplitude, there is only one weak phase in the 
SM, and therefore no direct CP asymmetry is expected. 
Any nonzero value from the experiment would be a signal 
for new physics. We define the relative size and strong 
phase of the two amplitudes as 



The BRs of B s — > Df^f^) can be estimated by using 
the experimental data on B s — > J/ipfoU^) together with 
the ratio of the BRs 



B(B° S -> D°f 



B(B° S 



Wo) 



Jo 



h 



B(B" -> J/^f 2 ) 



v cb v> 



Id 



VcbV* 



(1.3-1.5)%, 
0.8%, 



with the product of the form factors = (0.8 — 1.0) [12], 



Jo _ 



\V uh V: s A x /{V ch V: s A 2 )\ , 8 f ° = arg {A x /A 2 ) (10) 



Since both A\ and A 2 are from the same Feynman dia- 
grams, it is likely that Ai ~ A 2 , which implies ~ 0.5 
0. 

The neutral B s system is described by the mixing 



and 5g 



B a 

I2H ] and x f ^ s = 0.50 Q. A recent measurement [27[ 

B(B° S J/Mo) ~ B(B° S -+ 

- 0.2B(B° S -> J/%1>4>) - 2 x 10~ 4 

shows that the B s —> -D°/o(/2) decays have a BR of order 
10~ 6 . In this estimate the decays /o into ir + ir~ and f 2 
into K + K~ have been taken into account. 

It is straightforward to incorporate the B s de- 
cays into other light p-wave mesons, like /q(1370), 
fti(1170),/ii(1380), /i(1285) and /i(1420). But they re- 
quire high statistics to have an impact on 7, due to cither 
the suppressed production rates in B s decays [28| or the 
difficulty in the reconstruction of the decay modes [l(| ■ 

Finally, we remark on the BR estimate, which is ob- 
tained under the factorization approach in conjunction 
with the experimental data. The validity of this method 
can be tested by considering the ratios in the processes 
B° -> (D ,J/tP){K ,K* ): 



\B, 



-.p\B s )+ q \B s ),\BH)=p\B s )- q \B 



°) 



y K = B{B u -> D"K () )/B{B () 

with \p\ 2 + \q\ 2 — 1, and q/p denotes the weak phase in Vk* = B(B° — » D°K* )/B(B° 



the B s - B, 



ig q/p = V t * b V ts /(V tb V t * s ) 



> J/tpR°) ~ 1.4%, 
J/^K*°) ~ 0.5%,(13) 



In 



the SM, this ratio is close to unity and the phase (3 S is 
negligibly small (3 S ~ —0.019 rad. The normalized time- 
dependent decay widths are 0, [25| : 



r(B° s (t) D°(D°)f ) 



= p -t/r Bs 



1 



+ cos(Amt)C DO{DO)fo +sm(Amt)S DO{B o )fo ,(11) 

where T is the averaged decay width. For the correspond- 
ing B® decays, the plus signs in front of cosine and sine 
terms should be replaced by minus signs. Substituting 
the amplitudes defined in Eq. (J9j) , we have 

C D o fo = C DOfo = [1 - (r£) 2 ]/[l + (r f l)% 
S D o fo = -2r f l sin( 7 + 8* + 2(3 S )/[1 + (r f B °f], 
SDo fo = ~H° S M-l + S§. + 2A)/[1 + (r } B °f}. (12) 

The equality Cjjo f = Cqo f is a consequence of the 
uniqueness of the weak phase in decay amplitudes. Since 
both the strong phase difference 8g and the B s — B s 
mixing phase are expected small, Spo f Q and <5/jo j will 
have similar magnitudes but differ in sign. 



where the form factor products are used from Ref. [13l ]. 
Compared with the data tjk ~ 6.0% and yx* — 
3.2% 0, these ratios are theoretically undershot. If 
it is the same in B/B s decays into Kq 2 /(/oi f%)i the esti- 
mated BRs will be enhanced roughly by a factor of (4-6), 
which makes the proposed method more appealing. 

In summary, we have explored the possibility 
to extract the CP violation angle 7 with B — > 
(D°,D ,D° CP )K* (2) (U30) and B s -> (D°,D°)M (M = 
/o(980), / (1370), ^(1525), /i(1285), /i(1420), ^(1180)). 
A clean method is to use the two triangles formed by the 
decay amplitudes of B ± -> (D° ,D° ,D^ P )K*^ } (U30). 

We expect that B ± -> D ii'* ( ^ ) (1430) and 
5 ± — > D°K*^ (1430) have similar decay rates and 
the CP asymmetries have a strong correlation with 
7. Our method does not require the separation of 
the Cabibbo-suppressed D decays, which are usually 
buried under the combinatorial background. With the 
help of the factorization approach and the relevant 
experimental data we estimate the branching ratios of 
these modes to be of order 10~ 5 — 10~ 6 . Measurements 
of the branching fractions of Bd — > Di^^ 2 ^(1430) and 
time-dependent CP asymmetries in B s — > DM provide 
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an alternative way to extract the angle 7. No knowledge 
of the resonance structure in this method is required 
and therefore the angle 7 can be extracted without any 
hadronic uncertainty. 
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